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A  DYNAMIC  PROGRAMMING  SOLUTION  OP  A 
MISSILE  ALIDCATION  PROBLEM 

ABSTRACT 

A  generalization  of  a  missile  allocation  problem  proposed  by 
Piccariello^  is  formulated  using  dynamic  programming  techniques.  The 
formulation  leads  tc  an  efficient  algorithm  for  computing  integer  solutions 
to  the  discrete  problem. 
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INTRODUCTION 


H.  Piccarieilo  has  considered  the  following  interesting  missile  alloca¬ 
tion  problem.  Consider  a  target  complex  consisting  of  N  missile  launch 
siLes  and  M  control  centers.  It  is  assumed  that  if  one  (or  more)  control 
center  survives  an  attack,  then  all  surviving  missives  can  be  launched  and 
if  all  control  center.s  are  destroyed  chen  none  of  the  surviving  missiles  can 
be  launched.  Furthermore  the  assumption  is  made  that  each  of  the  N  +  M  missile 
sites  and  control  centers  presents  a  separate  target  to  an  enemy  missile 
attack,  -'his  is  equivalent  to  the  assumption  that  an  enemy  missile  can  destroy 
no  target  other  than  the  one  against  v;hich  it  is  directed.  Ihe  problem  is  to 
find  an  allocation  of  T  identical  missiles  that  when  attacking  the  target 
complex  will  .minimize  tne  expected  number  of  missiles  capable  of  being  launched. 

Piccarieilo  gives  a  solution  to  the  problem  when  it  is  considered  in  its 
continuous  form  (allocations  are  not  restricted  to  integer  values).  He  also 
investigates  the  discrete  form  of  the  problem  and  shows  by  example,  that,  in 
general,  the  solution  in  the  continuous  case  is  not  a  solution  for  the  discrete 
case.  He  then  derives  a  necessary  condition  for  the  existence  of  a  solution 
for  the  discrete  case  to  differ  from  the  continuous  case. 

In  this  paper  a  generalization  of  Piccarieilo 's  problem  is  formulated 
using  dynamic  programming  techniques.-  The  formulation  leads  to  an  efficient 
algorithm  for  computing  integer  solutions  to  the  discrete  problem. 

THE  FORMUIATION  OF  THE  PROBLEM 

We  shali  use  the  following  notations  and  definitions: 

probability  of  survival  of  missile  site  i;  i  =  1, N, 
when  attacked  by  a  single  incoming  missile, 

P.:.  probability  of  survival  of  control  center  j;  j  =  1,2,...M 
d  when  attacked  by  a  single  incoming  missile, 

a. :  number  of  missiles  at  missile  site  i, 

T;  total  number  of  attacking  missiles; 

E(X,Y):  expected  number  of  surviving  missiles  after  an  attack  by 

T  missiles,  where  Y  =  T  -  X  missiles  are  allocated  to  the 
control  centers  and  X  missiles  are  aliocu,ted  to  ‘•’he 
missile  sites. 


inie  following  expression  for  E(X,Y)  can  be  given;' 


E(X,Y) 


(1) 


where  X  is  the  number  of  missiles  allocated  to  the  ith  launcher  site,  is 
the  number  of  missiles  allocated  to  the  jth  control  center  and 

E 

)  X,  =X  (a) 


M 

Y’y.-Y^^'f  -  X.  (b) 

L. 

j=l 


The  objective'  is  to  find  the  minimum  values  of  E(X,Y)  subject  to  the 
constraints  a)  and  b)  where  each  x.  and  each  y^  is  an  integer. 

The  above  problem  reduces  to  the  Piccariello  problem  when  a^=a2=. . , 
Pl=Pg= . . .  =P..,  and  P^=P2= . . ..  =Pj,- 

FYNAMIC  PROGRAMMING  FORMULATION 

An  examination  of  equation  (l),  subject  to  the  constraints  (a)  and  (b) 
indicates  that  we  are  faced  with  an  N  +  M  dimensional  problem.  In  order  to 
circumvent  this  dimensional  difficulty.-  we  will  reformulate  the  problem  as  a 
multi-stage  problem  and  apply  the  techniques  of  dynamic  programming  to  obtain 

a  feasible  computational  scheme,  let  us  defin"", 

g  (x)  =  the  minimum  expected  number  of  surviving  missiles 
for  an  allocation  of  X  attacking  missiles  to  k  missile 
sites, 

and  ^  ^  minimum  value  of  the  probability  of  the  survival  of  at 

least  one  control  contcT  for  an  allocation  of  Y  a.-tacking  missxler 


to  i  control  centers.. 


Then  vze  have 


Min  E(X,Y)  =  Min  gjj(x)fj^(Y) 
{X,Y  j 

subject  to  the  constraint  X  +  Y  =  T,  where 

N  X. 

gjj(x)  =  Min  ^a.P.;  +...+  x^^  =  X, 

^^i ^  i=l 


(2) 


(5) 


ana 


1*1  V  . 

=  Min  1  -  TT  (1  -  p bubject  to  the  condition 


j-1 


'b'l 


1  2 

Using  t’ne  Principle  of  Optimality  we  can  express  the  g^j  recursively  as 


r  X 

gj^(x)  =  Min  I  aj^Pj^ 


,  k  =  2,3,...,N,  where  x^^  is  permitted 


to  vary  over  the  set  {0,1,2, .. .X).  For  k  -  1  we  have 


(5) 


gl(X)  -  a^P^ 


(6) 


To  obtain  fj^(Y)  we  first  observe  that  fj^(Y)  =  1  -  1^(y),  where  h  (y)  = 


M 


y . 

Max 


.Employing  the  Optimality  Principle  once  again,  we  obtain 


h^(Y)  =Maxr(l-P^^)  h^_^  (Y-yp  , 


-  2, 3 , . . . ,M  where  y .  e 

JCr 


{0,1,2, . . . ,Y}  and 


(7) 


Y 


(8) 


h^(Y)  -  (1  - 

The  value  Y  is  also  permitted  to  range  over  the  set  of  integers 

[0,1,2,. 


CCMPUTATIONAL  PROCEDURE 


The,  dynamic  programming  formulation  imbeds  the  original  problem  within 
a  family  of  analogous  problems  in  which  the  basic  parameters  N,  and  T 
assume  sets  of  values  which  permit  uj  to  obcaln,  in  the  course  of  the  compu¬ 
tation,  the  solution  to  a  variety  of  sub-problems.  Because  of  the  structure 
of  the  process  we  are  able  to  use  the  information  obtained  frcm  sub-problems 
of  the  originial  problem  to  obtain  th;  solution  of  the  original  problem.  The 
input  information  that  is  required  is  the  knowledge  of  K,  M,  and  T  and  the 
values  of  the  P.'s,  P.'s,  and  a.'s. 

'iv'e  begin  the  computation  by  obtaining  the  sequences  (gj^(X)}  and  {Xj^(x) ) 
from  equations  (5)  and  (6),  and  the  sequences  (h^(Y)J  and  (y^(Y) }  from 
equations  (7)  and  (8).  Given  this  information,  we  are  then  prepared  to  use 
equation  (2)  to  compute  Min  E(X, 


NUMERICAL  EXAMPLE 


Let  N  =  3,  M  =  2,  T  -  5,  P^^  -  .50,  P^  -  .30, 

a,  =  ]"  =  =  —  - 


and 


Pj_  =  .60,  Pp  = 

=  .50. 

g.(X)  =  a^P^, 

gp(x)  =  Min  « 

(v? 

g:.(x)  =  Min  - 
o  x^ 

h^(Y}  =  (1  - 

*1  )y 

h  (Y)  =  Max  • 

^  yo 

(U-. 

gp(X-x^)|-  ; 


=  .60,  a^ 

6,  7,  and 


=10 ,  ap  —  8 , 

8,  we  obtain  the 
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Using  the  previous  relations  we  obtain  the  values, 


g^(0)  =  10 

gl(l)  =  5 
£^(2)  =  2.5 
gl(3)  =  1-25 
g^(‘^)  =  .625 
g,(5)  =  .5125 


^(0)  =  a^  +  g^(0)  =  8  +  10  = 


/-»\ 

0(,  ..ou;  +  g^v-Ly  = 


18 


g^(lj  =  Min 


|_8(.50)^  + 


(0)  =  2.h  +  10  =  12.4 


g2(2)  =  Min 


8  +  g^(2)  =  10.5 
.4  =  gi(^)  =7.4 

.72  +  g^(0)  =  10.72 


=  7.4 


g.,(5)  =  Min  <v 


8  +  g^(5)  =  9-25 
2.4  +  g^(2)  =  4.9 
.72  +  g^(l)  =  5.72 

:.2l6  +  g^(0)  =  10.216  J 

8  +  g^(4)  =  8.625 
2.4  +  g^(3)  =  5.65 
(4)  =  Min  ^  .72  +  g^(2)  =  5.52 

.216  +  g^(l)  -  5.216 
.0648  +  gj^(O)  =  10.0648  J 


-  4.9 


=  5.52 


-  12.4 
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r 


62(5)  =  Min  J 


8  +  g^(5)  =  8.:j125 

2.4  +  g,(k)  =  5.025 

.72  +  gj^(3)  =  1.97 


g3(0)  = 


216  +  ^J2)  =  2.716 

0648  +  gj^(l)  =  5.0648 

^.01944  +  g^(o)  =  10.01944 

[12  +  30  4.  3]  ^50 


=  1.97 


12  +  g^il)  =  24.4 

g  (i)  =  Min  W 

5  ^  7.2  +  g^(0)  X  25.2  j 


{  =  24.4 


g,(2)  =  Min  < 


12  +  7.4  =  19.4 
7.2  +  12.4  =  19.6 
4.32  +  18  =  22.52 


r  =  19.4 


12  +  4.9  =  16.9 


g^(3)  =  Min 


7.2  +  7.i^  =  14.6 
4.32  +  12.4  =  16.72 
2.592  +  18  =  20.592 

12  +  5.52  =  15.52 


gj(4)  =  Min  4.32  +  7.4  ^  11.72 

2.592  ^  12.4  =  14.992 


>  =  14.6 


>  =  11.7.: 


^1.5552  +  18  =  19.5552^ 
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0 


=  Max  j 


> 


(.50)(,.8804)  =  .4402 
(.75)(.784)  =  .588 
(.875)(.64)  =  .56 
(.9375)(.40)  =  .575 
0 


We  are  finally  ready  to  make  use  of  eq.  (2) 

Min  E(Y,X)  =  Min  g„(x)  f  (Y) 

(X.Y)  ” 


=  .588 


where 


and 


1 


X  +  Y  =  T. 


For  oux  example  N  =  3,  M  =  2  so  that 


(X,Y} 

r 


Min  E(Y,X)  =  Min 
fX,Yj 


gjCx) 

fgCx) 

i.e.. 

fp{5) 

=  (50)(.4l2)  : 

=  12.56 

gjCi) 

=  (24.4){.52) 

=  12.688 

g:^(2.) 

1*2(5) 

=  (19.4)(.68) 

=  13.192 

gjC^) 

^2(2) 

=  (14.6)(.80) 

=  11.68 

s^C'O 

^2(1) 

=  (ll.Y2)(l)  = 

=  .11.72 

fo(0) 

C 

=  {9.22)(l)  = 

9.22 

=9.22 


The  minimum  value  is  9.22.  The  allocation  of  missiles  which  yields 
the  minimum  value  is  X  =  5,  Y  =  0.  What  remains  is  the  determination  of  the 


optimal  values  of  the  x^;  1  =  1,?,3: 

X.  =X^(X)  =x.(5)  =  y, 

Xg  =  x^(X-x_^)  =  x,^,(5)  =  1 

Tlitrefori-  thu  optiiQ-tl  allocation  poi  j  cy  is  p^iv*  n  by 
(x^,  x^.  x^,  y^,  y^)  =  (2,1, 2, 0,0). 


DIGCUSSION 


The  author  is  grateful  to  Mr.  ],'.  h.  Morrict  for  many  helpful  discussions 
during  the  fonaulation  and  solution  of  this  probiem.  Mr.  Mpi'ritt  has  pointed 
out  to  the  author  several  problem  areas  of  interest  in  which  this  procedure 
would  have,  direct  application.  For  in.stancc  it  would  be  possible  by  utilizing 
these  methods  to  decide  upon  an  optimiun  ratio  of  control  center  c,.  missile 
sites.  It  would  also  be  possible  to  use  this  methodology  to  lecermine  whether 
it  is  desirable  to  harden  ICBM  sites',  or  to  utilize  the  additional  money 
planned  to  be  spent  on  active  or  passive  defense  by  building  more  undefended 
launching  sites. 


la 


WIJJLIAM  SACC^ 
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